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We reveal that three-dimensional multi-orbital topological superconductivity can be identified by 
a bulk measurement, i.e., the temperature dependence of nuclear magnetic relaxation (NMR) rates. 

Below a critical temperature T c , the NMR rate in the topological state exhibits an anti-peak profile, 
which is opposite to the conventional s-wave state. This inversion coherence effect comes from a 
twist of order parameters with respect to orbital and spin degrees of freedom. Our self-consistent 
calculations in the model for CurcBhSea prove that the inverse coherence effect appears as a concave 
temperature dependence of the NMR rates. We propose that a time-reversal-invariant orbital- 
singlet spin-triplet topological superconductivity is characterized by the temperature dependence of 
the NMR rate. 

PACS numbers: 74.20.Rp, 74.25.Bt 


Physics and topology have mutual connections in mod¬ 
ern physics 0 Hj. Topological materials, such as topo¬ 
logical insulators and superconductors open an intrigu¬ 
ing avenue of materials science and quantum engineer¬ 
ing S0. A reliable way of detecting these materials 
prompts the study of topology in physics. Both surface 
and bulk probes are crucial for identifying topological 
materials. The bulk-boundary correspondence Q] indi¬ 
cates that gapless surface states between different topo¬ 
logical states signal nontrivial topological order. Bulk 
quantities also contain a signature of their topological 
states. A definite example is the conductance in integer 
quantum Hall systems; the quantized behaviors ruled by 
a topological invariant are observed Q. 


Topological superconductors have attracted intense at¬ 
tention. A classification study with different kinds of spa¬ 
tial dimension and symmetry Q indicates the presence 
of a topological superconductivity in three-dimensional 
(3D) time-reversal-invariant systems. Superconducting 
topological insulators, such as Cu^BUSea and 

Sni-zIn^Te [l2] are the candidates of the bulk 3D topo¬ 
logical superconductors. Their topology is argued by Z 2 
invariants Hill- A mean-field model of Cu^B^Sea [13] 
leads to notable predictions, such as vortex bound 
states [H, 0 and impurity effects 17 €3. Among these 
studies, although the gapless bound states on surfaces 
at low temperatures are actively studied [gl, [ll|, [2p| - l22| , 
observing an intrinsic behavior via bulk measurements is 
rarely argued. In addition, the setups of surface measure¬ 
ments would be sensitively affected by the characteristics 
of the interfaces between materials and probes. Thus, 
a bulk measurement sensitive to topological characters 
is important for supporting results in surface measure¬ 
ments and revealing the properties of Cooper-pair con¬ 
densations. 


A key quantity of connecting topological charac¬ 
ters with bulk measurements is a correlation function. 
Current-current correlation functions, for example, lead 


to the quantized conductance in integer quantum Hall ef¬ 
fects [5[. Spin-spin correlation functions are essential for 
superconductivity since these quantities well reflect the 
spin-state properties of Cooper pairs. Nuclear magnetic 
relaxation (NMR) rates (Tj _ ) are directly related to the 
spin-spin correlation functions. Therefore, it is interest¬ 
ing and important question to ask whether NMR rates 
contain any characteristic information on topological su¬ 
perconductivity. Specifically, the coherence effect is no¬ 
table. The NMR rate in the presence of a spin-singlet s- 
wave superconducting state is enhanced just below a crit¬ 
ical temperature T c , owing to the coherence factor p3 . 
This coherence peak (Hebel-Slichter peak) [H, [2f| comes 
from the formation of superconducting gaps. The ab¬ 
sence of the peak and a power-law behavior of T-j -1 at 
low temperatures indicate the occurrence of unconven¬ 
tional states nnn. 


In this paper we claim that the bulk measurements of 
NMR rates detect a 3D odd-parity fully-gapped topo¬ 
logical superconducting state in time-reversal-invariant 
multi-orbital systems. An inverse coherence effect just 
below T c is the signature of this odd-parity state; the 
coherence factor contributes to the NMR rates with 
an opposite sign to that of the conventional s-wave 
states. Our self-consistent calculations in the model of 
Cu^B^Sea [13 show that this sign reversal leads to an 
anti-peak behavior of the NMR rates below T c . Using 
the Fermion anticommutation property, we show that the 
odd parity allows a twist of a gap function with respect 
to the internal degrees of freedom (e.g., orbital and spin) 
in systems, causing the novel temperature dependence 
of NMR rates. Hence, NMR rates are bulk quantities 
sensitive to parity closely related to defining Z 2 invari¬ 
ants 0,113 in a superconducting state. 


We formulate the NMR rate of a multi-orbital su¬ 
perconductor, with a mean-field approach. The mean- 
field Bogoliubov-de Gennes (BdG) Hamiltonian is 7-L = 
(l/2)E k *Plm^ with i/’l = (4, c^) and 1 /£ = 





2 


(c£, c! fc ). The 2n 0 -component column (raw) vector Ck 
( 4 ) contains electron’s annihilation (creation) operators, 
with the number of orbitals n 0 . When n 0 = 2, we have 
c T = (cj, cj, c 2 , c 2 ). The BdG matrix is 




^o(fc) A(fe) \ 

At(fe) -m-k)}*) ■ 


(1) 


The normal-state Hamiltonian is Hq. The pairing po¬ 
tential (matrix) fulfills A T (— k) = —A(fc), owing to the 
Fermion anticommutation property. The check symbol 
(”) indicates a matrix in the Nambu space, whereas the 
hat symbol (~) does that in an orbital-spin space. The 
NMR rate [ME! 


is 


1 

Ti{T)T 



df(uj) 

duj 


x Re {p^'i^Pu'^) - p F ^'{u)[p F ™' M]*} . 

( 2 ) 


We use the unit system of H = fee = 1. The indices a and 
a' represent orbital labels. The Fermi-Dirac distribution 
function is denoted by f(u>). The spectral functions 
P G (oj) and p F (oj) are the submatrices of p G (oj) = 
(- 1 / 2 tt*) Efc [Gfc (zca n —t co -f- zO) — Gk(ico n —uj — zO)] • 
Temperature Green’s function is Gfc(zw ra ) = 
[iuj n — H’(fe)] -1 , with the fermionic Matsubara fre¬ 
quency uj n = nT(2n + 1) (n £ Z). The matrix form in 
the Nambu space is 


Gfe (icon') 


Gfc(icu n ) Ffc(zcu n ) 
Ffc(ica n ) Gfc(ica n ) 


( 3 ) 


The diagonal block, Gk leads to p G , relevant to electron’s 
density of states. The off-diagonal block, Fk contributes 
to the anomalous spectral function p F . 

The presence of a coherence peak just below T c is de¬ 
termined by the second term in Eq. © H3, Eli- The 
Hebel-Slichter peak appears when this term (including 
the minus sign in front of the spectral functions) has a 
positive contribution to Tj _1 - To intuitively understand 
the behaviors of the second term, we evaluate anomalous 
Green’s functions near T c . Linearizing G with respect to 
A, we obtain Ek 4(zw n ) « Efe G£ (iw n )A(k)G%(iw n ), 
with normal-state Green’s functions, GJ) and GJ). The 
spin-singlet property of s-wave states indicates A-jq = 
—Apr; we obtain —p F ^(p F ^)* = + \p F ^\ 2 , with p F ^ 0 . 
In contrast, the fc-integral with a spin-triplet p -wave gap 
leads to p F = 0 Q. 

Now, we apply the above arguments to topological 
superconducting insulators. We focus on the model of 
Cua;Bi 2 Se 3 . The normal electrons are effectively de¬ 
scribed by a two-orbital model (n 0 = 2 ), leading to a 


massive Dirac Hamiltonian, 


Ho(k) = 7° 


3 

pT° + ^ Vihr + mT 4 + h 5 (k)T 5 

i=1 


( 4 ) 


with chemical potential /lz, spin-orbit coupling constants 
Vi, and mass to. This Hamiltonian is the same as that 
in Ref. 03|, except / 15 . The last term corresponds to 
the effects of hexagonal warping in the Fermi surface 
of Cua,Bi 2 Se 3 [30]. We drop this term in most of our 
calculations, but argue the effects on the NMR rate at 
the end of this paper. Six kinds of 4 x 4 matrices T " 4 
(A = 0,1 ..., 5 ) ar e composed of the gamma matrices q^ 
(p = 0, 1, 2, 3)[3lj and the identity: = q A (A 7 = 4) 

and r 4 = I4, with q 5 = zq^'^q 3 - Our choice [32:]] is 
that 7 0 = a x <g> 1 2 , 7 1 = —icr y <g> s y , q 2 = ia y < 8 > s x , 
and q 3 = ia z <g> I 2 , where o x , y ,z {s x , y , z ) are the 2 x 2 
Pauli matrices in the orbital (spin) space. In this pa¬ 
per, we study the momentum-independent pair potential 
A, owing to the onsite interaction 0, [HI . The relation 
A T = —A leads to 


A = £ A^rVq 5 . 

A =0 


( 5 ) 


Since q 2 q 5 = I 2 < 8 > s y , the case of F " 4 to be the iden¬ 
tity (i.e., A = 4) describes a spin-singlet s-wave state: 
A 4 oc E( c -fc|4t + c -fci4t>- The additional (A ^ 4) 
characterizes a twist of each order-parameter component 
in the orbital-spin space, compared to the conventional 
s-wave state. According to Ref. 03}, even-parity or¬ 
der parameters (Ai g states) are given by A 4 and A°. 
Odd-parity states correspond to A 1,2,3,5 . The compo¬ 
nent A 5 corresponds to an ocld-parity fully-gapped (Ai„) 
state 0, E3: A 5 oc E(c_ fci 4 t + c 2 _ k ^c l ki ). The odd- 
parity state on Dirac-type normal electrons given by 
Eq. © 33 Ms nontrivial, in terms of not only a Z 2 
invariant 1, l]|, but also a ID Z 2 invariant 0, 34]. 

Let us study the anomalous spectral function p F of the 
odd-parity state near T c . Normal-state Green’s functions 
are evaluated by an algebraic relation of Hq; we find that 
[H' 0 (k)] 2 = [E{k)] 2 , with H' 0 = H 0 + p. This property 
corresponds to the fact that the Dirac equation is the 
square root of the Klein-Gordon eauation[3l}. Hence, 
we obtain G k (iu) n ) = Ez=± Pz{k)/[iuj n - £E(k) + p\, 
with P± = (1/2) [1 ± (Hq/E)]. The projectors P( are 
rewritten by Pe(k) = 7 ° Ea=o w t(k)T A . Then, when 
A (k) = A 5 = A 5 F 5 q 2 q 5 and /15 = 0 , we find that 


P 4; A 5 ) = x (w)A 5 + % (w)q u A 5 , 


( 6 ) 


near T c . The coefficients are X 5 ’ 40 = 

Efc T,e l' Ne,e'(k,uj)Wg^°{k), where Nt,t> = 
{l/[{t + l')E - 2 p}[6(u> - IE + p) - 6(u + PE - p)], 

= ELo 4 »E 44 and W 4 °< = w$w%-w° t w%. 
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FIG. 1. (Color online) Temperature dependence of nu¬ 
clear magnetic relaxation rates in (a) an odd-parity gap 
A 5 oc XXc-fcjCfc f + c 2 _ hr ci j.) ( A lu ;i3 or A 2 [ft] ) and (b) 
an even-parity gap A 4 oc ]T( c -fc4. c fet + c -fet c fct) 0 or 

Ai a @] )• We set the chemical potential p = 0.8 and the 
Dirac mass m = 0.4 [See, Eq. ©]. We tune either U or V so 
that the zero-temperature gap amplitude takes 0.01. Inset: 
temperature dependence of the superconducting gaps. 


Since T 5 7 2 7 5 = —ICJy ® A 5 is odd under an orbital- 
index exchange, while is even under a spin-index 
exchange. The multiplication of j° with A 5 does not 
change the property of a spin-index exchange. Thus, 
we show that p F “ a ' (w) = p F “ a '(u:). This property 
indicates that the second term in Eq. © negatively 
contributes to T-j -1 : — p F ^ a ‘ (p F ^ a ')* = — \p F ^ a '\ 2 . An 
inverse coherence effect can occur in the presence of the 
odd-parity state. 

Now, we show the inverse coherence effect of T-j -1 , solv¬ 
ing the BdG equations self-consistently. We focus on 
Eq. © with /15 = 0. We set Vi = v, for simplicity. 
We take a unit system with v = 1. We vary to, with 
fixed chemical potential p = 0.8. The superconduct¬ 
ing pairing model is two-orbital on-site density-density 
interaction IH, Hiil; #mt(*) = U [ni(x ) 2 + n2(x) 2 ] + 
2Vni(x)n2(x) with n a = i c s ^ c s- The gap equa¬ 

tion is 

Aaa' = — T aa ' n aa i, (7) 


with h F = — Ylk lim T ->o-i- exp(— ivj n T)F k {iuj n ) and a = 
(a, s). The intra-orbital coupling constant U and the 
inter-orbital coupling constant V are bundled up by V aQ /, 
i.e., V aa ' = U with a = a' and V aQ / = V with a ^ a'. 


Self-consistently solving Eq. 0, we evaluate the temper¬ 
ature dependence of Eq. 0 - The k integrals are per¬ 
formed by the trapezoidal rule in the spherical coordi¬ 
nate system, with cutoff momentum fc max = 9 and mesh 
(iVfc, Ng, Nf) = (128,128,128). The smearing factor of 
the delta function is set by 0.01. We tune either U or V 
so that the zero-temperature gap amplitude takes 0 . 01 . 
In this paper, we focus on the odd-parity gap function 
A 5 oc S( c -fc4, c fct + c -fct c fc.|,) an< ^ t ' ie even-parity gap 
function A 4 oc ^(c 1 i, | c L- + c 2 _ k ^c k ^). Assessing the lin¬ 
earized gap equation [13| , we find that the former appears 
as the most stable phase when U = 0 and V < 0 (attrac¬ 
tive inter-orbital interaction). The latter occurs when 
U < 0 and V = 0 (attractive intra-orbital interaction). 





FIG. 2. (Color online) Nuclear magnetic relaxation rates of 
an odd-parity gap A 5 ( A\ u [h| or A 2 3 ), changing the 
Dirac mass m with fixed chemical potential p = 0.8. The 
total rate (a), the normal part (b), and the anomalous part 
(c) are shown. 

Figure |T| shows the temperature dependence of 
(TiT) -1 , with to = 0.4. When the gap function is 
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odd parity [Fig. Q3a)], an anti-peak or concave behav¬ 
ior occurs below T c , coming from the anomalous part of 
(T^T) -1 (blue triangles). The temperature dependence 
is contrast to that of Fig. |TJb) , in which the gap function 
is even parity and the coherence peak appears just below 
T c . Thus, we obtain the inverse coherence effect in A 5 . It 
is worth noting that the temperature dependence of the 
gap function (insets in Fig. [T|) has no significant differ¬ 
ence between A 5 and A 4 . Moreover, we can find that the 
exponential temperature dependence of {T\T)~ l occurs 
at lower temperatures; this is consistent with the fully- 
gaped feature of A 5 . Thus, the temperature dependence 
of the NMR rate definitely signals the emergence of A 5 . 



m/p 


FIG. 3. (Color online) Height of coherence peaks [normal¬ 
ized by the value at T c , Ti(T c )T c ], changing m/p (Dirac 
mass/chemical potential). The negative height means the in¬ 
verse coherence effects. 


Let us examine how the inverse coherence effect de¬ 
pends on material variables, focusing on mass m. Fig¬ 
ure m a) shows the m-dependence of ( T\T ) 1 when the 
gap function is A 5 . The concave behavior is much pro¬ 
nounced when m is small. The normal parts [Fig. [2jb)] 
of the NMR rate do not considerably depend on m, while 
the anomalous one [Fig. Oc)] goes to zero when m —> p. 
Let us evaluate the extreme value of the anomalous part 
below T c normalized by Ti(T c )T c . This quantity char¬ 
acterizes the amount of the coherence effect. Figure [3] 
shows the extrema (i.e., height), changing m/p. We find 
that when we take A 4 , the extrema is positive and al¬ 
most independent of m/p. In contrast, in the case of 
A 5 the extrema is negative and drastically depends on 
m/p. The value approaches zero when m/p —> 1, while it 
takes a negative constant when m/p —>• 0. When m ~ p, 
the odd-parity fully-gaped gap function is mapped to the 
spin-triplet p-wave gap function with d-vector d oc fc[ 3. 
Thus, the inverse coherence effects are suppressed when 
m becomes large. We also find that the negative extrema 
as m/p —> 1 is proportional to an indicator of relativis¬ 
tic effects [18]], /3 = |ufcF|/|m| = 1 / (p/m) 2 — 1, where fcp 
is the Fermi momentum. It indicates that a relativistic 
effect, which is intrinsic in the Dirac Hamiltonian, char¬ 
acterizes the inverse coherence effect of the NMR rate. 
We stress that the anti-peak behaviors are predominate 


when m/p —> 0. According to Refs. 17 -T 9 }, the odd- 
parity state with small Dirac mass is robust against non¬ 
magnetic impurities. Thus, one can detect the inverse 
coherence effect even in dirty samples. 


Now, we propose an unified point of view on NMR 
rates near T c , depending on different gap functions. 
We focus on three kinds of parity p s , p m and p 0 on 
A. Spin parity p s is related to the exchange of spin: 
A“,“ (fc) = p s A““ (fc) with p s = ±1. Similarly, mo¬ 
mentum (orbital) parity p m (p 0 ) is determined by ex¬ 
change of momentum (orbital): A“^f (—fc) = p m A“^f (fc) 
[ A“^, a (fc) = PoA““’(fc) ] with p m(o) = ±1. The fun¬ 
damental relation of A is A T (— fc) = —A (fc), leading to 
PsPmPo = —1- They are not independent of parity; in 
the present model the odd orbital parity is forbidden in 
the even-parity fully-gapped states (H} . Table I summa¬ 
rizes (p s , Pm, Po) of different gaps. The non-zero contri¬ 
butions of p F to Tj _1 is allowed when p m =l. In fast, a 
single-orbital spin-triplet p-wave state (p m = — 1) has no 
peak [U 0]. The sign of the coherence effect is deter¬ 
mined by the spin parity p s . A multi-orbital odd-parity 
state allows both p m = 1 and p s = 1 [Tij]. 

We argue the application of NMR rates to finding 3D 
topological superconductors. The knowledge of normal- 
state Fermi surfaces is needed [I, 35j, to pick up plausible 
samples in terms of Zi invariants. After such a screening 
test, NMR measurements would lead to one of the evi¬ 
dence for the odd-parity fully-gapped state. The details 
are summarized in Ref. 


Finally, we discuss the NMR rate with anisotropic odd- 
parity states labeled by E u : A 1 oc W (c 1 ,^c}^ + c 2 | cj | ) 
lldl ll9L 130} . According to Ref. [30|, this pairing can be 
induced in the presence of a hexagonal warping term 
/15 = *A(fc3. + fcl) with k± = k x ± iky. Near T c , we 
can find that p F ^ aa (w; A 1 ) oc A. Since the orbital parity 
is odd but the spin parity is even, the inverse coherence 
effect can occur within |A| ^ 0. We will conclusively 
discuss the coherence effect of the NMR rate in topologi¬ 
cal superconductors elsewhere, including anisotropy and 
momentum dependence in the gap functions. 


In summary, we showed a bulk measurement, the 
NMR rate near T c as a sign of topological supercon¬ 
ductivity. The inverse coherence effect was predicted 
in the NMR rate below T c . Our self-consistent calcu¬ 
lations in the model of CmcB^Sea lead to the temper¬ 
ature dependence of (TjT) , with a concave behav¬ 
ior below T c . This notable effect originates from the 
anomalous part of the NMR rate in the odd-parity gap 
A 5 oc Y^( c2 -ki c k-\ + c2 -k\ c k])- Moreover, the negative 
height when m/p —> 1 is proportional to an indica¬ 
tor of normal-electron relativistic effects. Therefore, we 
claim that the detection of the inverse coherence effect 
definitely signals the emergence of bulk topological su¬ 
perconductors. Our approach of focusing on correla¬ 
tion functions unveiled a link of bulk measurements with 
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TABLE I. Parity table of different gap functions. The coherence effect is characterized by spin parity p B [A““ (k) = p E (fc)], 
momentum parity p m [A“|* (fe) = p m A““ (—fe)], and orbital parity p a [A““ (fc) = p 0 A“j“(fe)], with p B p m p 0 = —1. 


Gap type 

Spin parity p s 

Momentum parity p m 

Orbital parity p 0 

Coherence effect 

2-orbital s-wave SC. (Ai 9 ) 

-1 

+ 1 

+ 1 

positive 

2-orbital topo. SC. (Ai„) 

+ 1 

+ 1 

-1 

negative 

1-orbital s-wave SC. 

-1 

+ 1 

+ 1 

positive 

1-orbital p x + ip y -wave SC. 

+ 1 

-1 

+ 1 

0 


a topological superconducting state. Accordingly, the 
present study can be a clue of building an explicit formula 
of topological invariants, with aid of bulk-measurement 
quantities in superconductors. 
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onal warping term. The calculations were performed by 
the supercomputing system PRIMERGY BX900 at the 
Japan Atomic Energy Agency. This study was partially 
supported by JSPS KAKENHI Grant Number 26800197 
and Grant-in-Aid for Scientific Research (S) Grant Num¬ 
ber 23226019. 
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Supplemental material 

We show a way of classifying time-reversal-invariant (TRI) 3D topological superconductors, in terms of three 
kinds of parity. The setting of physical model here is more general than that in the main text. We impose that 
normal electrons are subjected to a two-orbital Hamiltonian with time-reversal and spatial-inversion symmetries, 
and the inversion transformation, P is expressed by one of the orbital Puali matrices, say, a x . Three kinds of 
parity are related to, respectively, orbital degrees, spatial inversion, and Fermi-surface topology [l|. The orbital 
parity of pairing potential A(fc) is defined by partial transpose with respect to orbital indices, given by 7^[A(fc)] = 
(l/2)[A(fc) + J2i=x y,z( a i ® l)A(fe)(crj (g) 1)*]. Thus, we have 70[A(fc)] = p 0 A(k), with p 0 = ±1. A straightforward 
calculation shows that this definition is equal to that in the main text. The spatial-inversion parity of A(fc) is defined 
by P t A(fe)P* = pA(-fc), with p = ±1 ill]. The parity of Fermi-surface topology is given by [1| 

pf =n sgn t e26 ( ri )] sgn [ £ 2fc( r j)]> (!) 

b 

where C 2 b(Fj) is the eigenvalue of normal-electron Hamiltonian Hq^T^ at TRI momentum point Fj in the first Brillouin 
zone. Since the system has time-reversal symmetry, the Kramers degeneracy at indicates that £ 2 b(Fj) = £2&+i(Fj)- 
Thus, the subscript b runs over half of labels in the eigenvalues of Ho{k). According to Ref. [H, within weak-coupling 
superconductivity, a ID Z 2 invariant of a TRI odd-parity state is defined by (—lj'd'W'l =_Dp, with a closed path CV, 
connecting between and Tj. If (—l) ! d c *d = _l ; the system is topologically nontrivial [l|. Moreover, odd numbers 
of zero-energy states are predicted on the boundaries made of opening up a closed path CV,. Table [TT] shows all the 
eight patterns of (p 0 , P, Pf) and the corresponding ID Z 2 when p = —1. We show that, in a two-orbital system, cases 
V and VI are forbidden since p 0 and p are not independent of each other 0 • 

Now, we summarize a method of detecting TRI 3D odd-parity topological superconductors, using the link of the 
classification table with the inverse coherence effects in nuclear magnetic relaxation (NMR) rates. If the gap function 
is an odd function of momentum (e.g., oc k x + ik y ), no coherence peak is expectable. Therefore, we focus on the 
case of even momentum parity in Table I of the main text; the inverse coherence peaks can appear when p Q = — 1 . 
Cases VII and VM in table ITT1 allow an anti-peak behavior in NMR rates. We stress that, without empolying any NMR 
experiments in a superconducting state, case VH is distinguished from case VM by the Fermi-surface topology attainable 
from normal-state probes, such as the first-principles calculations of energy band and angular-resolved photo-emission 
spectroscopy of Fermi surfaces. After a screening test of Fermi surface topology, NMR rate measurements should be 
performed. Therefore, the observation of inverse coherence peaks is the evidence of the 3D odd-parity fully-gapped 
topological superconductivity. 


TABLE II. Parity classification of time-reversal-invariant 3D superconductors. 


Case 

Orbital 

Spatial inversion 

Fermi-surface topology 

ID Z 2 

I 

+ 1 

+ 1 

+ 1 

- 

II 

+ 1 

+ 1 

-1 

- 

n 

+ 1 

-1 

+ 1 

even 

IV 

+ 1 

-1 

-1 

odd 

V (forbidden) 

-1 

+ 1 

+ 1 

- 

VI (forbidden) 

-1 

+ 1 

-1 

- 

vn 

-1 

-1 

+ 1 

even 

vm 

-1 

-1 

-1 

odd 


[1] M. Sato, Phys. Rev. B 81, 220504(R) (2010). 

[2] Using the linear map of partial transpose with respect to orbital degrees, we find for each A A (oc r A 7 1 2 * 7 5 ) that 2 p a = 

1 + p — m y + m z , with m y , z = ±1. This identity indicates that odd orbital parity is not compatible with even parity. 







